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The infinite color-Coulomb energy of color-charged states is related to enhanced density of near- 
zero modes of the Faddeev-Popov operator in Coulomb gauge. We confirm the enhancement 
in numerical simulations and show that it is tied to the presence of percolating center vortex 
configurations. 
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1. Introduction: Confinement scenario in Coulomb gauge 

The Faddeev-Popov operator in Coulomb gauge 

M(A) = -V-^(A), where ^ C (A) = dj8 ac + f abc A b i (x), (1.1) 

plays a crucial role in the Gribov-Zwanziger confinement scenario [[j], ||]. It enters the kernel K in 
the (classical) Coulomb energy of a color charge distribution p : 

H coul = \J d 3 xd 3 y p"(x)K« b (x,y;A)p b (y), p a = pLuer ~ ■ E c , (1.2) 

where 

K ab (x,y;A) = [M~ l (-V 2 )^ 1 ]^. (1.3) 

The essence of the scenario can be summarized in the following way: The Coulomb-gauge 
condition V ■ A" = does not fix the gauge completely. Gribov [[[]] suggested to restrict to the 
subspace of transverse gauge fields for which the Faddeev-Popov operator is positive, i.e. local 
minima with respect to g(x) of 

I[A,g}= Jdx[^A a (x)]\ where 8 A i = g- l A i g + g- l d i g. (1.4) 

The boundary of this Gribov region (GR) is called the Gribov horizon. However, even this does not 
eliminate the Coulomb-gauge ambiguities completely, one has to further narrow the gauge-field 
configuration space to the fundamental modular region (FMR), i.e. the set of absolute minima of 



the functional (1.4). Both the GR and the FMR are bounded in every direction and convex. The 
dimension of the gauge-field configuration space is huge, so it is reasonable to expect that most 
configurations are located close to its boundary (horizon). The interaction kernel K contains the 
inverse of the FP operator, which is strictly zero on the horizon and near-zero close to the horizon. 
A high density of configurations near the horizon leads to a strong enhancement of the Coulomb 
interaction energy, and hopefully causes color confinement. 

In this contribution we formulate a simple criterion of confinement for static color charges 
through properties of eigenstates of the FP operator in Coulomb gauge close to the Gribov horizon, 
and then discuss how the fulfillment of this criterion depends on presence/absence of center vor- 
tices. Details, as well as some analytic insights on the connections between center vortices and the 
Gribov horizon, can be found in a recent publication [gj]. 1 

2. Lattice Faddeev-Popov operator and its eigenstates 

If we parametrize link variables in SU(2) lattice gauge theory by 

U li {x) = b il {x) + ia c a c il {x), b^x) 2 + £a c ^x) 2 = 1, (2.1) 



'We have also investigated localization properties of the lowest nontrivial eigenvectors of the Faddeev-Popov oper- 
ator in Coulomb gauge. These were discussed in Jeff Greensite's talk at this conference and in Sect. V of Ref. H. 
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the lattice Faddeev-Popov operator in Coulomb gauge is given by the following expression: 

M% = 8« b £U n ,[b k ( x )+b k ( x -k)]-8 xv _Mx)-^-My)} 
k ^ } 

- ^LfVKW - VKW} • (2-2) 

k k J 
We are interested in its eigenstates 

L< = K 0i" )fl , (2.3) 

in particular in their properties near to the Gribov horizon (i.e. in the limit X — » 0). The most 
relevant quantities are 

• the density of eigenstates p{X), and 

• the average Laplacian: 

f^l0i' l)a (-V 2 )xv^" )a *. (2-4) 

a,xy 

3. A confinement condition 

We shall now formulate a simple confinement criterion in terms of properties of eigenstates of 
the FP operator. The energy of a static color charge state *P"[A;x] in Coulomb gauge 

g = %*^* c) ~ WolHcoM ~ (K""(x,x;A)) (3.1) 

can be easily shown to be given by 

1 — . / F \ r*™* dX 

*=W^W g0ing, ° L X^<^» fo<V^~. (3.2, 

An immediate consequence is that the excitation energy <§ of a static, unscreened color charge 
is divergent if, at infinite volume, 

ltm mpm >0 . 03) 

A^O X 

This criterion is a necessary but not sufficient condition for confinement; an explicit example will 
be given at the end of Sect. ||| (It is obviously not fulfilled in the free theory, where p(A) ~ V~X, 
F(X) = X, and consequently S' ~ \/X max .) 

4. Three ensembles of lattice configurations 

We will investigate fulfillment of the condition ( |3~3| ) in three ensembles of configurations: 

1. full configurations, {^(x)}; 

2. "vortex-only" configurations: these are obtained from full configurations fixed to the (di- 
rect) maximal center gauge [||] by center projection, {Z M (x) = sign 

3. "vortex-removed" configurations, obtained by the recipe of de Forcrand and D'Elia [^]: 

{U^\x)=zl{x)U li {x)}. 

Each configuration in these three ensembles was brought to Coulomb gauge by maximizing 
with respect to gauge transformations, on each time slice & cou i(t) = Y*xY,l=i jTr[Uk(x,t)]. 
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Figure 1: p(A) and F(X) for full lattice configurations. 



5. Results 

Pure gauge theory at zero temperature. The results for full configurations at j8 = 2.1 are 
shown in Figure [I], for a series of lattice volumes. 2 Both p(A) and F(A) exhibit a sharp "bend" 
near A — > 0, and behave near like a small power of A . A scaling analysis similar to that used in 
random matrix theory gives the estimates 

p(A)~A a25 , F(A)~A 038 . (5.1) 

The confinement condition ( |3.3[ ) is obviously satisfied, which is a direct manifestation of the mech- 
anism proposed by Gribov and Zwanziger. 

The situation in vortex-only configurations is displayed in Figure ||. The enhancement of the 
density of states is even more pronounced than in full configurations, and both quantities of interest 
seem to converge to a non-zero value in the infinite volume limit 

p(0)~0.06, F(0)~1.0. (5.2) 

(though their proportionality to very small powers of A cannot be excluded). Once again, the 
condition (|0|) is fulfilled. 



Eigenvalue density, p=2.1 



F(J.)=<<p l I -V 2 I <p x >, p=2.1 
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Figure 2: p(A) and F(X) for vortex-only (center-projected) configurations. 



2 The results for jS = 2.3 and 2.4 can be found in Ref. [H], and are qualitatively the same as those for j8 =2.1. 
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Eigenvalue density, p=2.1 F(X) s ( <p x -V 2 I <p x >, |3=2.1 
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Figure 3: p(X) andF(A) for vortex-removed configurations. 

The eigenvalue spectrum of the FP operator is drastically different for vortex-removed config- 
urations, see Figure || for the largest available, 20 4 lattice. The density displays a series of peaks, 
and values of F(X) are organized into bands, separated by gaps. This can be understood rather 
simply: For the Laplacian operator (equal to the FP operator at zero-th order in the gauge coupling) 
the eigenvalue density, at finite volume, is a sum of delta-functions, and each eigenvalue is multi- 
ply degenerate. The vortex-removed configuration seems to be just a small perturbation around the 
zero-coupling limit, which lifts the degeneracy. In this way, delta-functions in the density of states 
turn into distinct peaks of finite width, and degenerate values of F(X) spread into bands. The num- 
ber of values inside the k-th band of F(X) exactly matches the degeneracy of the k-th eigenvalue of 
the unperturbed Laplacian operator. 

This result demonstrates a deep relation between the Gribov-horizon and center-vortex con- 
finement mechanism. Center vortices seem to be the field configurations providing the mechanism 
needed for enhancement of eigenvalues near the horizon. 

Pure gauge theory in the deconfined phase. A seemingly paradoxical result is obtained above 
the deconfinement transition: our quantities above T c look the same as at T = (cf. Figure || with 
Figs. H] and ^ ! However, one should keep in mind that spacelike links are a confining ensemble 
even in the deconfined phase, and spacelike Wilson loops have an area law behavior. 

The result for the deconfined phase can be naturally explained in the Gribov-horizon scenario. 



Eigenvalue density, p=2.3, T=2 F(X) = < <p x I -V 2 | <p x >, p=2.3, T=2 




Figure 4: p (X ) and F(X) in the deconfined phase, for full and vortex-removed configurations. 
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In Coulomb gauge the gauge fixing is done independently on each 3d time slice. According to the 
horizon scenario, on each time slice, 3d configurations A(x) are favored that lie near the horizon 
of a 3d gauge theory, and this enhances the instantaneous color-Coulomb potential. This is true for 
every temperature T, including in the deconfined phase, because temperature determines the extent 
of the lattice in the fourth dimension. Thus, the horizon scenario provides a framework in which 
confinement may be understood, but it is not detailed enough to tell us under what conditions the 
infinite color-Coulomb potential may be screened to give a finite self-energy. 3 

6. Conclusions 

The low-lying eigenvalues of the FP operator in Coulomb gauge tend towards zero as the 
lattice volume increases. The density of the eigenvalues goes as a small power of A, and this, 
together with a similar behavior of the average Laplacian, F(X), assures the infrared divergence 
of the energy of an unscreened color charge. These facts support the ideas of the Gribov-horizon 
confinement scenario. 

The constant density of low-lying eigenvalues can be attributed to the vortex component of 
gauge-field configurations. A thermalized configuration in a pure gauge theory factors into a con- 
fining piece (the vortex-only part), and a piece which closely resembles the lattice of a gauge-Higgs 
theory in the Higgs phase (the vortex-removed configuration). This establishes firm connection be- 
tween the center-vortex picture and the Gribov-horizon scenario. 

The Gribov-Zwanziger scenario, though invented to explain confinement, is operative also in 
the finite temperature deconfined phase. 

Here we only covered results of our numerical investigations. Related analytical developments 
were omitted and can be found in our recent publication [^J|. 
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